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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
research MEMORANDUM 

PRELIMINARY STUDY OF SOME FACTORS WHICH AFFECT THE 
STALL-FLUTTER CHARACTERISTICS OF THIN WINGS 
By A. Gerald Rainey 

SUMMARY 


The results of an exploratory, analytical, and experimental study 
of some of the factors which might he of importance in the stall flutter 
of thin wings are presented. The factors considered were Mach number, 
Reynolds number, density, aspect ratio, sweepback, structural damping, 
location of torsion nodal line, and concentrated tip weights. 

The importance of aerodynamic torsional damping in the stall flutter 
of thin wings has been demonstrated through comparison of regions of 
negative torsional damping and regions of flutter. 

The results of a series of experiments on a thin wing tested at 
various lengths Indicate that compressibility definitely alters the 
stall-flutter characteristics of wings of this type and that the com- 
pressibility effects appear to depend on the aspect ratio. 

A brief study of the Inertia effects of concentrated tip weights 
Indicates that such effects can be important; however, the large number 
of parameters involved makes it difficult to generalize the results. 

An approximate analysis is presented for obtaining an estimate of the 
stall-flutter characteristics of particular wing -weight configurations. 

Some of the other factors considered were found to be more or less 
significant; however, all the factors would require further study if 
their effects are to be more fiilly interpreted. 


INTRODUCTION 


In general, as the angle of attack of a wing is increased to values 
near the angle of stall, the flutter velocity is decreased to values 
much lower than that observed at 'angles near zero lift. Since the 
potential-flow theories of classical flutter do hot Indicate an effect 
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of lift coefficient, the decrease in flutter velocity is presumed to be 
associated with nonpotential flow and has come to be known as stall 
flutter even though the wing may not be completely stalled in the usual 
sense of the word. Wear zero lift the type of flutter which occurs is 
usually a coupled, bending-torsion type, whereas flutter encountered at 
the higher angles of attack is usually predominantly a single-degree- 
of -freedom flutter occurring in the torsional mode. 

Much information exists concerning- the stall flutter of moderately 
thick wings at low speeds (refs. 1 to 6) and the effects of several 
parameters on the stall flutter of thin propeller blades have been 
presented in reference 7- While stall flutter has not been considered 
a serious problem in the design of conventional wings (previous to 1950), 
experience with thin propeller blades indicates that stall flutter may 
become an important consideration in the design of thin wings. Since 
designers are considering the use of thin wings for high-speed airplanes, 
it was considered desirable to reexamine the stall-flutter problem. 

The purpose of this paper is to present the results of an explora- 
tory, analytical, and experimental study of some of the factors which 
might be of importance in the stall flutter of thin wings. The factors 
considered are Mach n\unber, Reynolds number, density, aspect ratio, 
sweepback, structural damping, location of torsion nodal line, and 
concentrated tip weights. The results obtained by variation of any 
particular parameter are somewhat sketchy since the purpose of this 
investigation is to search for the significant parameters rather than 
to define completely the effect of the significant parameters. 


SYMBOIS 


a 

9-0 

9-n 

b 

Ca 

c 

dt 


speed of sound, ft/sec 

initial amplitude of oscillation, radians 
amplitude of oscillation after n cycles, radians 
semichord perpendicular to leading edge, ft 
torsion spring constant per unit span, ft-lb/ft-radian 
chord, ft (2b) 

time derivative of energy, ft-lb/sec 
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fa(x) 

^‘hi 

fai 

f 

exp 

g 

Sa 

la 

Ka 

L 

M 

Ma 

% 


n 


( 


rof- 


K 


R 

S 

V 

X 


spanwise torsional deflection function 
first natural bending frequency, cps 

first natural torsion frequency, cps 

experimental flutter frequency, cps 
damping coefficient 

structural damping coefficient for torsion mode 

moment of inertia per unit span about the elastic axis, 
sl\ag-ft^/ ft 

spring constant of damping apparatus, ft-lb/radian 
wing length, ft 

twisting moment per unit span, ft-lb/ft 

complex aerodynamic moment per unit span per radian deflec- 
tion ft-lb/ft-radian 

t 

aerodynamic moment per radian in phase with velocity, / 

ft-lb/radian 

coefficient of aerodynamic moment per radian in phase with 
velocity 

number of cycles 

effective inertia parameter where r^^ is the effective radius 
of gyration and k is the effective mass density ratio 

Reynolds mmiber . 

wing area, ft2 

velocity, ft/sec 

distance along span, ft ■ 

angular deflection, radians 


a 
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angular velocity, radians/sec 

angular acceleration, radians/sec^ 

angle of attack at root, degrees 

amplitude of torsional oscillation, radians 

angle of sweepback, degrees 

density of test medium, slugs/ft^ 

frequency of oscillation, radians/sec 

first natural torsional frequency, radians/sec = 2nfa^'' 

METHODS AND APPARATUS 


• The experimental investigation consisted of essentially two phases: 
(l) conventional flutter tests with the addition of the angle of attack 
as a primary variable and (2) aerodynamic damping measiirements with the 
angle of attack again being a primary variable. 


Wind Tunnel 

The experiments were conducted in the Langley 4.5-foot flutter 
research tunnel which is of the closed-throat single-return type 
employing both air and Freon 12 as a testing medium at pressures from 
one atmosphere down to l/8 atmosphere. During the course of this 
investigation the test section of the 4.5-foot-diameter flutter tunnel 
was altered to a rectangular section 2 feet by 4 feet to accommodate 
other types of investigations. The effects of this change in test 
section have not been evaluated but they are believed to be small. 


Models 

Five thin wings were used for the experiments and are illustrated 
in figure 1. Wing number 1 was originally 24 inches longj however, in 
order to investigate the effects of aspect ratio, this wing was shortened 
in steps of three inches and the various configurations are designated 
wing number 1-24, 1-21, l-l8, and so forth. When equipped with a tip 
weight the wing designation is followed by the letter W, such as 1-21W. 
The model designations, properties, and test conditions are listed in 
table I. 
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Wing munber 1 vas of solid aliiminum-alloy construction and was 
used in studying the effects of Mach number^ aspect ratio, and con- 
centrated weights on stall flutter. This wing also was used to measure 
the aerodynamic damping moment for torsional oscillations. A series of 
foirr wings (designated 2, 3^ and 5) varying in sweep angle and • 
aspect ratio were employed to study the stall-flutter characteristics 
of swept wings. The construction of these wings consisted of wood 
laminations on a stainless steel plate and the length-to-chord ratio 
was kept constant at 2.7- 

Each of the five wings was equipped with resistance-wire strain 
gages arranged to obtain a record of the flutter frequency and to 
indicate the relative magnitude of the strains in the bending and 
torsion degrees of freedom at flutter. 


Flutter Tests 

During the flutter tests the models were mounted as cantilever 
beams on a heavy, rigid raoimt which could be rotated and clamped at 
any desired angle of attack. With the wing clamped at the desired 
angle of attack, the airspeed of the' tiinnel was slowly increased until 
a steady or divergent oscillation of the model occurred. At this point 
the tunnel conditions were observed and an oscillographic record of 
the strain-gage outputs obtained. The velocity at this point was con- 
sidered to be the minimum flutter velocity for the wing tested at that 
particular angle of attack. In some cases the amplitudes encountered 
at this velocity were sufficiently mild so that the velocity could be 
increased further until some limiting condition was reached, such as 
excessive vibration amplitudes, excessive static bending load, or 
maximum tunnel velocity. The same procedure was then repeated for 
angle -of -attack increments of 2^ or less at various angles throughout 
the range from 0^ to 2^+-^. 

The effects of density, Reynolds number, and Mach number were 
studied by performing the above experiments at various pressures from 
atmospheric down to l/8 atmosphere in both air and Freon 12. 

The effects of aspect ratio and Mach number were investigated by 
shortening wing number 1 in steps of three inches. 

The effects of concentrated tip weights were studied, briefly, by 
clamping a solid ”bendalloy" ellipsoid 10 inches long and 0.75 inch 
in diameter to the tip of wings number 1-21 and l-l8. The center of 
gravity of the weight was placed at the 50 -percent- chord line, 0.375 inch 
inboard from the tip in the plane of the wing. The characteristics of 
the wings equipped with the weight are listed- in table I under the 
designation 1-21W and l-l8W. 
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Aerodynamic Damping Measurements 
6 

The aerodynamic torsional damping moment was measured by the method 
of decaying oscillations (ref. 8). 

A special wing mount (shown schematically in fig. 2) was designed 
and built for wing number 1-24 to allow torsional oscillations but 
restrict the bending mode. The root of the wing was restrained in 
torsion only by two relatively soft coil springs. The first natural 
torsional frequency of the wing-shaft-spring system was about l/8 of the 
first natural torsional frequency of wing number 1-24 when mounted as a 
cantilever. In the special mount the wing tip was restrained in bending 
by a steel cable and, consequently, the first natural bending frequency 
of the wing in the damping mount was appreciably higher than in the 
cantilever condition. For these reasons the wing was considered as a 
rigid body (constant angle of attack spanwise) when oscillating at the 
first natural torsional frequency of the wing-shaft-spring system. 

The damping of the system was obtained by deflecting the wing- 
shaft-spring system in torsion to an amplitude of approximately 3° from 
the mean angle of attack. The wing was then released and the resulting 
oscillation was recorded. The damping coefficient was computed from 
the following relationship (ref. 9): 


_ 1 ®'0 

® Jtn a^^ 


and the aerodynamic damping moment was computed from 


^ (^total " Stare) ^ “ Saero^ 

where gtotal damping coefficient obtained with air flow and 

Stare the damping coefficient obtained with no air flow. The non- 
dimensional aerodynamic damping moment coefficient is defined as 


m. = — q-..-. 

^V^Sc 

2 

The aerodynamic torsional- damping moment coefficient was obtained 
at various angles of attack by rotation of the entire wing-shaft-spring 
system to the desired angle of attack. The coefficient was determined 
as a function of v/bo) for angles of attack of 0°, 4°, 8°, 12°, l6°, 
and 20°. " " 
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The effects of the location of the axis of rotation or torsion 
nodal line were studied by clamping the wing to the shaft so that the 
axis of the shaft coincided with the desired chordwise location. The 
nodal line positions studied were 25 , 32 . 5 , 50, and 75 percent of the 
chord rearward of the leading edge. 

The maximum Mach number reached in the damping measurements was 
about 0.2) the maximum Reynolds number was about 0.8 X 10^. 


RESULTS AKD DISCUSSION ■ 

General Remarks Concerning the Effects of Lift 
Coefficient on the Flutter of Wings 


In general, it has been found that as the angle of maximum lift is 
approached the, flutter velocity is reduced appreciably from the value 
obtained at angles near zero lift. In addition, the mode of flutter at 
the higher angles of *attack is usually predominantly a single—degree— 
of -freedom flutter occurring in the torsional mode. Since this type of 
^"l^"tter is associated with nonpotential flow^ the phenomenon, as pre- 
viously noted, has come to be known as stall flutter. 

A typical low Mach nimiber variation of the flutter velocity coef- 
ficient with angle of attack for thin wings is shown in figure 3« The 
flutter velocity is plotted in the nondlmensional form V/bo^ S'S a 
function of the angle of attack at the root for wing number 2. The 
lower boundary in figure 3 represents the velocity at which an oscilla- 
tion first started. The upper curve from = 6° to = l8° 
represents the velocity at which. the oscillation reached dangerous 
amplitudes. The area labeled "flutter free" between 20° and 2h° is 
bounded by the velocity coefficients at which the oscillation stopped. 

In this range of angle of attack the velocity was increased until the 
static bending load on the wing became excessive. The limiting velocity 
is indicated by the upper dashed curve. 


As the angle of attack is increased the flutter velocity coefficient 
decreases rapidly and reaches a minimum value near = l6°. The 

^®lo®fty coefficient v/bo^ at this minimum point is approxi- 
mately equal to 1.0. The velocities listed in table I correspond to 
this minimum value for the various configurations and test conditions 
shown. Experimental evidence obtained on numerous propeller blades in 
reference 7 and on several thin wings in the present investigation 
indicates that the rule-of-thumb, V/bcwQ^ — 1.0, for estimating the 
minimum flutter velocity of simple thin wings may be quite useful. Of 
course, the exceptions to this elementary nile are numerous as will be 
i^^i®9-t6d in the subsequent discuss ion o f the various factors studied. 
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The remainder of the data shown in figure 3, other than the minimum 
point, do not lend themselves to generalities other than, perhaps, a 
few speculative remarks regarding the possible causes of the flutter. 

At zero angle of attack the flutter was a coupled bending-torsion type 
and the flutter frequency was at some value between the natiiral first 
bending frequency and the natiiral first torsional frequency. . As the 
angle of attack was increased, the flutter— velocity coefficient decreased 
rapidly and the flutter frequency gradually approached the torsional 
frequency. This decrease in flutter velocity is probably due to a 
decrease in aerodynamic damping resulting from partial separation and 
boundary -layer time-lag effects. McCullough and Gault in reference 10 
indicate that, for steady flow, relatively sharp-edged airfoils exhibit 
a meas\rrable region of separated flow near the leading edge at angles 
of attack as low as a = 1°. As the angle of attack is increased this 
separated-flow region is extended and, presumably, the boundary-layer 
time-lag effects become more pronounced. The mechanism of this decrease 
in damping is not clearly understood- and would be difficult to treat 
analytically. 

As- the angle of attack is increased to values well beyond the 
steady-state stall angle (ai = 20 ° to. 2 l)-°), the vibration encountered 
resembles a forced vibration rather than a self-excited vibration. The 
amplitude of vibration increased gradually as the velocity was increased 
beyond the initial vibration point then, after reaching a maximum value, 
decreased to zero at the velocities Indicated in figure 3. In this 
region of angle- of attack, the vortex frequency, estimated on the basis 
of the projected wing area, was of the same order as the vibration 
frequency which was within a few percent of the natural torsional 
frequency. Although the motion under these conditions closely resembles 
a forced oscillation, it is difficult to separate the phenomenon at very 
high angles from the phenomenon at moderate angles and for that reason 
the oscillations at all angles of attack will be referred to as flutter. 


Some Considerations of Aerodynamic Damping 

As long ago as 1928 Frazer and Duncan (ref. l) hypothesized that 
the decrease in flutter speed with increasing angle of attack may be 
accounted for by a decrease of the aerodynamic torsional damping. In 
1943 Mary Victory (ref. 5 ) used a conventional flutter analysis altering 
only the torsional damping coefficients to calculate with good agree- 
ment several stall-flutter conditions obtained experimentally by Studer 
(ref. 2). The damping coefficients used in Victory's analysis were 
obtained experimentally on a similar wing at about the same values of 
Reynolds number and frequency as the flutter data of Studer. 

Some of the aerodynamic torsional damping coefficients measured 
using wing number 1-24 are shown in figure 4(a). Contours of constant 
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values of the coefficient are shovm in figure Ma) as functions of 
V/bo) and angle of attack with the aixis of rotation at the midchord. 

The curves above the zero damping line represent increasing values of 
negative damping, except at the higher angles of attack where the damping 
decreased to zero then became positive again at higher values of V/boj. 
The area enclosed by the zero damping curves represents a region where 
a single-degree-of-freedom torsional instability is possible. 

For comparison, the actual flutter regions obtained for the same 
wing (no. 1-2^4-) while mounted as a cantilever are shown in figure 4(b). 
One interesting feature of these data lies in the similarity of the 
flutter regions of figure 4(b) to the regions of negative damping shown 
in figvire 4(a). This similarity is rather striking when it is realized 
that the damping data were obtained when the wing was moiinted in a 
mechanism which provided appreciably different boundary conditions from 
those provided by the simple cantilever mount. The frequency of the 
wing moimted in the damping mechanism was about one-eighth of the 
torsional frequency of the wing mo\mted as a cantilever. The damping 
mechanism severely restricted the bending degree of freedom^ whereas 
the cantilever mount, of coirrse, allowed this degree of freedom. The 
maximum Mach number of the damping measurements was about 0.2 as com- 
pared with a Mach number of about 0.75 for the flutter tests. The 
dissimilarity between the two sets of data at the lower angles of 
attack is attributed to compressibility effects which will be- discussed 
later. In view of the similarity of the regions which occurred in spite 
of the widely different conditions which applied to the two determina- 
tions, it appears that the torsional aerodynamic damping is very impor- 
tant in the stall flutter of thin wings. 


Stall-Flutter Analysis 

The problem of calcixlatlng stall-flutter velocities has been a 
subject of much interest and some of the more promising approaches are 
discussed in reference 6. At the present time, no satisfactory stall- 
flutter analysis has been developed. This lack of an adequate analytical 
approach is due, primarily, to the difficulty in treating the air forces. 
The problem is somewhat similar to that of predicting the- maximum lift 
coefficient. Presumably, if' the oscillatory air forces were known for 
the proper range of significant variables, the analysis of stall flutter 
would be relatively simple. Unfortunately, the large number of signifi- 
cant variables makes it impractical to obtain adequate air-force coef- 
ficients except over very limited ranges. 

Observations of the stall-flutter characteristics of the unswept 
wings used in the present investigation indicate the possibility of 
applying a relatively simple analysis employing the measirred aerodynamic 
torsional damping coefficients. The observations are: (l) The flutter 
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occurs in almost a pure torsional mode with very little translation of 
the elastic axis and (2) the flutter frequency is very nearly equal to 
the first natural torsional frequency. (See table I.) Based on these 
observations a simplified single-degree-of-freedom analysis (developed 
in the appendix) yields the following relation: 



V Y 

n \bcoJ 


= 0 


Since the structural damping is always positive, at flutter the. aero- 
dynamic torsional damping must be negative and at least equal to the 
structural torsional damping. 


If the aerodynamic torsional damping moment coefficients are 
available for the proper range of significant variables, the flutter 
velocity can be calculated. The aerodynamic damping moment coef- 
ficients measured in the present investigation are presented in fig- 
tire 5 for the various angles of attack and axes of rotation. ’ The 


(the form most con- 


2nw/ y Y 

figure shows the values of the quantity 

n \^J 

venlent for the calculations) as functions of v/bcn. For a given wing 
the quantity 



may be estimated, calculated, or measured, and the flutter .velocity 
may be determined from the figures at each angle of attack where the 
aerodynamic damping is negative. Of coiirse, at the lower angles of 
attack where the flutter can no longer be represented as a single- 
degree-of-freedom torsional vibration, the analysis is meaningless. 

' Effects of structural damping .- One of the parameters which can 
cause deviations from the "rule of thumb" mentioned previously, that 
is, ( V/b<%.)mj^jj^ = 1.0, is the structural damping coefficient g^^. This 

effect of struct\iral damping (based on the preceding analysis) is 
illustrated in figure 6 where the minimum flutter velocity is shown 
as a function of the effective inertia parameter for several values 
of structural damping. The analytical results shown are for the axis 
of rotation at the midchord. The horizontal line for zero structural 
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damping indicates that the minim\am flutter -velocity coefficient wovild 
be independent of the inertia parameter for this hypothetical condition. 
The addition of a small amount of damping increases the minimum flutter 
velocity coefficient appreciably and makes the coefficient a function 


of 



A thin solid aluminum-alloy wing at sea-level density 


would have a value of 



near 50. Increasing the structviral 


damping coefficient from O.OQl to 0.010 would approximately double the 
minimum flutter velocity coefficient. 


Effects of fluid density . - Variation of the density can cause 
deviations from the nominal minimum flutter velocity coefficient of 1.0. 
This. effect is illustrated in figure 7 where the minimum flutter velo- 


city coefficient is plotted as a fimction of the quantity 

which is inversely proportional to density. The solid curve, which 
was calculated for the midchord axis of rotation by the analysis 

(r 2\ 

previously discussed, rises sharply at values of gaj-^) representing 

high density then tends to level off at' values corresponding tb reduced 
density. The hypothetical thin aluminum-alloy wing mentioned earlier 

would have a value of of about 0.05 for a damping coefficient 


of Sa 0.001 at sea level. The analysis indicates that the minimum 


flutter velocity coefficient would be increased from about 1.0 at sea 
level to 1.5 by decreasing the density to that at 40,000 feet altitude. 


Effects of location of the axis of rotation .- The effects of the 
location of the axis of rotation can be estimated with the aid of the 
aerodynamic torsional damping coefficients for various axes of rotation. 
Examination of figure 5 indicates a definite reduction in the extent of 
the negative damping regions for the forward locations of the axis of 
rotation. It appears, therefore, that the stall-flutter problem for 
wings having a torsion nodal line well forward may be greatly relieved 
as compared with the problem for wings having rearward locations of the 
torsion nodal line. 


Effects of Reynolds Number 

The effects of Reynolds number on stall flutter have not been 
thoroughly investigated; however, a few speculative remarks are possible. 
At low Mach numbers the aerodynamic characteristics of thin unswept wings 
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at high angles of attack are relatively insensitive to changes in 
Reynolds number. For this reason, it seems that Reynolds ntunber might 
be relatively ^important in the stall flutter of thin vinswept wings 
at low Mach numbers. This supposition is a least partially verified 
by some of the data shown in figure 7 where it can be seen that the 
minimum flutter-speed coefficient for wing number 1-24 increased from 

_X_ = 0.750 to = 0.846 over a range of density. The corre- 

bo^ bo^ 

spending range of Reynolds number was from 3*4l X 10° to O.085 x 10°. 
The increase in flutter velocity coefficient with decreased density is 
in qualitative agreement with the trend predicted by the analysis. 

The effects of Reynolds number at high Mach mxmbers and for swept 
wings are appreciably more obscure and require further investigation. 


Effects of Mach Number and Aspect Ratio on the Stall Flutter 

of Unswept Wings 

The "rule -of -thumb" mentioned earlier, that is, = I’O 

\ ^/min 

or (V)ni -tn = could be used to predict the minimum Mach number at 

which flutter might occur if compressibility did not alter this minimum 
flutter velocity coefficient. The minimum flutter Mach number 
would simply be b<%/a> however. Baker (ref. 7) found that, for thin 
propeller blades, (v/bo:a)inin ^ function of Mach number. By varying 

the speed of sound a, it was shown that the minimum flutter Mach number 
became a nonlinear fvinction of bcqx/a at the higher Mach ntunbers, and 
at sufficiently high values of bixct/a (approx. 0.5) no flutter occurred 
up to the limit of the tests (l^ip “ 1-4). 

The results of the present investigation of compress lb ility effects 
are shown in figure 8. The minimirra flutter Mach number is shown as a 
function of the quantity bc^/S' wings number 1-24 to 1-15 • Since 

it was not possible to obtain a sufficiently wide variation of bcq^/a 
by variations of the speed of sound, it was necessary to increase (%, 
by reducing the length of the original wing number 1-24. Consequently, 
the effects obtained may be thought of as combined effects of Mach 
niunber and aspect ratio. The results obtained on wing number 1-24 show 
a fairly linear relationship between the minimum flutter Mach number 
and bo^a and there was no indication of an approaching limiting value 
of boqx/a beyond which flutter would not occur. The value of boq^/a 

reached in these tests was about 0.43 which is slightly less than the 
value of 0.5 presented in reference 7 as the value required to eliminate 
stall flutter of propeller blades. Wing number 1-21 was tested at values 
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of ' bc%/a- up to 0.5 and^ although flutter was not eliminated at any 
value of ha^/ a^ the shape of the curve of minimiim flutter Mach number 
as a function of ba^/a (shown in fig. 8) indicated that a limiting 

bo^ 

value was approached at = 0.5* With wing number I-I 8 a new curve 

a 

bo^ 

was defined with an apparent limiting value somewhere near = 0.59* 

a 

bo^. 

When wing number 1-15 vas tested at =‘0.7, no flutter was obtained 

a ’ 

up to the choking Mach number of the tunnel (M = 0.8); however^ inter- 

bcjo 

mittent flutter did occur at =. 0.67 indicating that a more com- 

a 

plete curve for this wing might have a shape similar to the dashed 
curve shown. 

The results of this series of experiments on one thin wing at 
various lengths indicate that compressibility definitely alters the 
stall-flutter characteristics of wings of this type. Furthermore^ the 
compressibility effects appear to depend^ to a large extent^ on* the 
aspect ratio. 


Effects of Sweepback 

The available data on the effects of sweepback on stall flutter of 
thin wings are limited to a. series of experiments on wings number 2 ^ 3 ^ 
and 5 - In figure 9 the variation of the flutter velocity coefficient 
with angle of attack at atmospheric pressure is presented for the four 
wings varying in angle of sweepback from 0 ° to 45° • 

The curves for the 0^^ 15°^ and 30° wings are very similar indi- 
cating no significant effect' of sweep angle up to 30 ^ for this series 
of wings. The 45^ ving, however^ did not experience flutter at high 
angles of attack even though the velocity was increased to about 80 per- 
cent of the flutter velocity obtained at zero angle of attack. The 
tests on the 45° wing were stopped at conditions indicated by the dashed 
curve because of the excessive static bending load on the wing. 

The results of the present experiments indicate the possible 
importance of sweepback; however^ additional research Will be required 
before the sweepback effects can be applied to other configurations and 
test conditions. 
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Inertia Effects of Concentrated Tip Weights 

The inertia effects of concentrated tip weights on the stall- 
flutter characteristics of thin wings were investigated by equipping 
wings number 1-21 and l-l8 with a dense ellipsoid. High-density 
material ( "bendalloy") was used to obtain a small voliune so that the 
aerodynamic effects of the mass would be a minimum. The aerodynamic 
effects of large bodies on stall flutter have not been investigated. 

The weight reduced the natviral torsional frequency of both wings approxi- 
mately by a factor of 0.5- 

The variation of the flutter velocity coefficient with angle of 
attack for wings number 1-21 and l-l8 both with and without the weight 
is shown in figure 10. When tested without the weight both wings had 
a minimvim flutter velocity coefficient of about 1.0, and this minimum 
occurred between = l6° and = 20°. When tested with the weight, 

the minimum flutter velocity coefficient was increased to about 3.0 and 
the angle of attack at which this minimum occurred was decreased to 
about = 8° or 10°. For comparison, the calculated variation (see 

section entitled "Stall-Flutter Analysis"), for both wings with and 
without the weight are shown in the same figure. The calciilated results 
show the same trend as the experimental results in that the calculations 
indicate a large increase in the minimum flutter velocity coefficient 
.and a decrease of the angle of attack at which this minimum occurs when 
the wing is equipped with a weight. 

The results shown in figure 10 indicate that the minimum flutter 
velocity coefficient for these two wings was approximately tripled by 
equipping the wings with a concentrated tip weight, and the minimum 
flutter velocity was approximately doubled. It should not be concluded, 
however, that the addition of a weight always increases the minimum 
flutter velocity. The analysis on which the calculated results are 
based indicates that the minimum flutter velocity for a thin unswept 
wing at low Mach numbers depends on the amount and distribution of mass 
or inertia, the fluid density, and the structural damping. For example, 
when wing number l-l8W was tested in Freon (high density) the minimum 
flutter velocity coefficient ' was about 1.5 as compared with 1.0 for wing 
number l-l8. The actual velocity, however, was a few percent lower for 
the weighted condition. 


Because of the many variables of Importance in the stall flutter 
of wings with concentrated tip weights^ it is not possible to generalize 
the results of this brief investigation. The inertia effects of a 
particular wing -weight configuration can be estimated by use of figure 7 
where the minimum flutter velocity coefficient is shown as a function 


of the quantity 



The symbols represent experimentally 
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determined values for several wing-weight-density combinations and the 
solid curve represents calculated minimum values for wings having a 
torsion nodal line near the midchord. The analysis agrees qualitatively 
with the experiments. 


CONCLUDING REMARKS 


The importance of aerodynamic torsional damping in the stall 
flutter of thin wings has been demonstrated through comparison of 
regions of negative torsional damping and" regions of flutter. 

The results of a series of experiments on a thin wing tested at 
various lengths indicate that compressibility definitely alters the 
stall-flutter characteristics of wings of this type and that the com- 
pressibility effects appear to depend on the aspect ratio. 

A brief study of the inertia effects of concentrated tip weights 
indicates that such effects can be important; however, the large nxnnber 
of parameters involved makes it difficult to generalize the results. 

An approximate analysis is presented for obtaining an estimate of the 
stall-flutter characteristics of particular wing-weight configurations. 

Some of the other factors considered (Reynolds number, sweepback, 
fluid density, structural damping, and location of torsion nodal line) 
were fovmd to be more or less significant; however, all the factors 
would require fixrthsr study if their effects are to be more fully 
interpreted . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. 
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APPENDIX 


. A simplified single-degree-of -freedom-flutter analysis using 
experimentally determined aerodynamic damping coefficients has been 
developed in order to provide some means for estimation of stall- 
flutter velocities of thin wings. . • 

The basic assvunptions of the analysis are as follows: 

1. The motion of the wing is restricted to a single degree of 
freedom, namely, torsion. 

2. The flutter frequency cjo is assumed to be equal to the natirral 
torsional frequency ob,. 

If the approach used in reference 11 is followed, the total moment 
contributed by an element of length dx is 


M = Ia,a + Ca(l + iga)a + MaO. (l) 


vrhere Ma is the complex aerodynamic moment per luilt span per radian 
deflection. - 

The time rate of change of energy of each of the elements is equal 
to its total moment multiplied by its angular velocity. The rate of 
change of total energy for the entire wing having only torsional motion 
is simply the sum of the energy rates contributed by each of the 
elements 


( 2 ) 


dt 


Md dx 


At the flutter condition the time rate of change of energy must equal 
zero. Substituting the expression for the moment, equation (l), into 
the expression for energy, equation (2),, and setting the total rate of 
change of energy equal to zero produces 


[^d + Ca(l + iga)cL + Mg^c^d. dx = 0 


dE 

dt 


'L 


(3) 
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The motion, of each element is harmonic and can be expressed as 


a 


= 



eioit 


where represents the spanwise mode shape for torsional vibra- 

tions. It is assumed that the aerodynamic terms do not appreciably 
change the mode shape. 

If this relation is used and Ca = is defined^ equation ( 3 ) 

becomes 








la |?a(x)l ^dx + 



(4) 


Equation (4) can be separated into a real and an imaginary equation 
both of which must vanish to provide a^ solution to (4). Since it has 
been observed that in most cases of stall flutter cd = cOa^ the real 
part of equation (4) must approach zero. Normally the real part of 
equation (4) is used to provide the flutter frequency but, since the 
flutter frequency has been ass\nned to be equal to. the natioral torsional 
frequency, the solution of the real part of equation (4) will not be 
required. 


The imaginary part of equation (4), solutions of which provide the 
flutter velocity, is 




Ia[fa(x3 ^dx + 


(imag. )Ma [fa(x]] ^dx = 0 (5) 


or 


%.^ga Ia[fa(x)J^dx + 2pV^ md b^ [fa( xj] ^dx = 0 (6) 


The first part of equation (6) represents the structural damping 
energy whereas the second part represents the aerodynamic damping energy. 
'If it is assumed that a representative value of the semichord can be 
chosen for purposes of defining v/b(%^ that is, ir^ not a function 
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of X, then equation (6) can be simplified as 


Ia[£a(x)]^dx 


. = 0 


n Iboi 


V 


K /e ^ It \^ci^ 


where represents the effective or mean value of the important 

r^2 

inertia parameter — 


C 
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Figxare 1.- Flutter models. 
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(a) Axis of rotation at 0.25c. 

2mA /Y V 

Figure 5.- Variation of — ^ — ) with — for various angles of attack. 
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(b) Axis of rotation at 0.325c. 


Figure Continued. 
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(c) Axis of rotation at 0.50c. 


Figure 5.- Continued. 
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Figure 7.- Comparison of experimental and calculated minimum flutter 
velocity coefficients for wings nxmiber 1-24 to 1-18 at various 
densities with and without a concentrated tip wei^t. 
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No flutter up to choking Mach number 
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Angle of attack, ac^ 


Figure 10.- Comparison of experimental and calculated flutter velocity- 
coefficients for wings number 1-21 and 1-18 with and without concen- 
trated tip wei^t at l/8 atmosphere in air. 




